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ABSTRACT 

The  main  objective  of  the  present  study  is  to  provide  experimental  information 
about  turbulence  structures  in  unsteady  flow  fields. 

Measurement s  were  made  in  the  wake  of  a  NACA  0012  airfoil  oscillating  in  pitch 
around  an  axis  at  25%  chord  distance  from  the  leading  edge  with  a  sinusoidal 
motion.  Instantaneous  profiles  of  velocity  and  all  non-zero  Reynolds  stresses  are 
determined  at  48  phase  angles  within  the  oscillation  period.  Measurements  were 
made  at  .7,  2#  5,  10  and  20  %  chord  downstream  of  the  trailing  edge.  Various 

incidences>  from  zero  to  just  above  stall  limit*  various  amplitudes  and  frequencies 
were  considered.  The  reduced  frequency  (^c/2U)  ranges  from  0  to  1.  2  at  a  chord 
Reynolds  number  of  300*000. 

An  improved  rotating  slanted  hot  wire  anemometer  technique  is  used  for  the 
determination  of  all  non-zero  Reynolds  stresses. 

A  micro  processor  controlled  system  samples*  digitizes  and  records  on  tape  the  hot 
wire  anemometer  signal  as  well  as  the  wing  and  probe  position  as  functions  of  time. 
Instantaneous  profiles  are  reconstructed  through  ensemble  and  conditional  averaging 
techniques. 

Results  are  presented  and  compared  with  steady  state  values.  Detailed  results  are 
made  available  on  magnetic  tape. 
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d i  incident*  increase 

E  hot  wire  anemometer  voltage 

H  wake  shape  factor 

i  instantaneous  incidence 

k  chord  frequency  coefficient  <^c/2G„> 

m  Bruun's  power 

Q  absolute  velocity 

Qeff  effective  cooling  velocity 

r  distance  from  oscillation  axis 

s  Reynolds  stress 

u  meanst ream  velocity  fluctuation 

v  velocity  fluctuation  normal  to  meanstream  and  span 

w  spanwise  velocity  fluctuation 

Vc  chordwise  velocity 

Vp  pitchwise  velocity 

y  pitchwise  coordinate 

aQ  wire  slanting  angle 

aP*  probe  rotation  angle 

6  wake  displacement  thickness 

e  wake  momentum  thickness 

l  probe  inclination  angle 

0  flow  angle#  phase  angle 

0  yaw  angle 

^  rms  of  fluctuations 

(periodic)  time  averaged  value 
(r)  relative 


1  HffROPUCHON 


Unsteady  phenomena  appear  in  many  f luid-dynamic  systems  as  a  basic  event  or  'as  a  _ 
perturbation  of  a  steady  state  situation.  In  external  aerodynamics  for  example# 
the  flow  along  airfoils  is  often  strongly  unsteady  due  to  the  motion  or  the 
vibrations  (flutter)  of  the  blades  (e. g.  helicopter  blades#  propeller  flows# 
Darr ieus-type  windmills).  In  such  situations  the  increasing  demand  for  higher 
performance  requires  a  more  basic  understanding  of  the  unsteady  phenomena  in  order 
to  be  able  to  improve  the  existing  design  methods. 

Although  many  aspects  of  unsteady  flows  can  be  handled  by  linear#  or  non-linear 
inviscid  theory#  important  classes  of  phenomena  are  dominated  by#  or  strongly 
dependent  on  viscous  effects  such  as  dynamic  stall#  stall  flutter  or  the  turbulence 
structure  of  wakes  from  rotating  blades  in  turbomach ines#  which  can  strongly 
influence  the  flow  on  the  downstream  stator. 

The  actual  spectacular  development  of  computational  methods  to  solve  steady  and 
unsteady  flows  has  reached  a  stage  where  two-dimensional#  time  averaged  Naviei — 
Stokes  equations  can  be  solved  within  reasonable  cost  and  storage  requirements#  for 
both  steady  and  unsteady  situations  C303.  The  comparison  of  calculated  Naviei — 
Stokes  solutions  and  experimental  data  is  however  still  limited  mainly  by  the  lack 
of  detailed  knowledge  of  the  turbulence  input  to  the  equations  C201.  More 
particularly  in  the  case  of  unsteady  boundary  layers#  the  question  of  the  influence 
of  the  unsteadiness  on  the  turbulence  structure  is  still  unanswered. 

At  the  present  time  several  experiments  are  undertaken  in  this  field  :  experi¬ 

ments  on  flat  plates  C6-9# 14# 15#  23-26#  313#  airfoils  C17#27,  283  or  dlffutors  C293 
where  the  boundary  layers  are  exposed  to  an  oscillating  pratsure  gradient.  Several 
experiments  are  made  on  unsteady  jets  C323  and  on  pitching  or  plunging  airfoils 
C4# 12# 13# 10# 19# 221.  In  the  latter  case  experiments  are  mainly  dealing  with 
unsteady  airloads  C4#221#  flow  visual isation  C18# 191  and  to  some  extend  with 
velocity  end  turbulence  distributions  C12# 131. 
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The  objective  of  our  research  is  therefore  to  provide  detailed  information  about 
the  turbulence  structure  in  makes  and  boundary  layers  of  an  oscillating  airfoil  by 
comparing  mean  velocities  and  all  components  of  the  Reynolds  stress  tensor  with  the 
corresponding  steady  state  values 

In  order  to  concentrate  on  the  turbulence  properties*  the  research  mill  be  limited 
to  incompressible  flouts  without  heat  transfer. 

In  the  present  report*  the  measurement  and  data  reduction  techniques  are  detailed 
and  all  available  results  of  the  experimental  work  done  in  the  wake  of  the 
considered  airfoil  are  presented  and  discussed. 


2  PREVIOUS  UORK 


Some  excellent  reviews  on  the  subject  of  unsteady  flows  can  be  mentioned  the 
Freeman  lecture  by  lie  Croskey  C213  which  is  devoted  to  unsteady  flows  in  general* 
reviews  of  Telionis  C333  which  are  mainly  devoted  to  unsteady  boundary  layers  and  a 
recent  review  of  experimental  work  done  in  this  field  by  L.  W.  Carr  E33.  The 
authors  insist  on  the  urgent  need  for  experimental  data  on  turbulence  effects  in 
unsteady  flows  both  as  basic  data  -  comparison  with  calculations  and  turbulence 
models  as  well  as  providing  insight  into  the  physical  mechanism  of  energy  exchange 
between  the  random  turbulent  motion  and  the  organized  periodic  flow. 

In  the  following*  more  details  are  given  about  past  work  done  which  is  closely 
related  to  the  present  experiments*  namely  work  done  on  unsteady  boundary  layers 
and/or  wakes  on  flat  plates  and  airfoils. 

Most  of  the  flat  plate  experiments  are  done  on  a  fixed  plate  with  an  external 
oscillating  flow  Detailed  experiments  are  available  with  and  without  adverse 
pressure  gradient  in  unseparated  and  to  some  extend  in  separated  flows.  Frequency 
coefficients  based  on  the  boundary  layer  thickness  range  from  O  to  .24  1233  and 
oscillation  amplitudes  reach  507.  of  the  outer  flow  161.  Data  of  velocities*  two  or 
three  Reynolds  stresses*  static  pressure  and  skin  friction  are  available  Table  1 
gives  an  overview  of  work  done  in  this  field  with  some  details  about  each 
■e-xp-ejilment.  This  table  is  not  believed  to  be  complete  but  it  gives  an  idea  of  the 
type  of  e  x  pVf  rmvnt-s  ^eod  results  available  in  the  open  literature. 


The  flat  plate  experiments  can  be  considered  as  a  more  fundamental  research  whereas 
more  direct  data  are  needed  about  current  flow  types  such  as  flows  around  common 
used  airfoils  Experiments  on  unsteady  flows  around  airfoils  are  in  general 
performed  using  an  oscillating  airfoil  (pitching  or  plunging)  in  a  steady  free 
stream*  whereas  the  flat  plates  are  generally  fixed.  As  experiments  with  a  fixed 
airfoil  in  an  oscillating  external  flow  one  can  mention  the  work  of  Satyarayana 
C273  and  of  Saxena  et  al  C283.  Saxena  et  al  made  experiments  at  angles  of  attack 
close  to  the  stall  limit.  In  these  experiments  pressure  distributions  and  velocity 
profiles  are  presented*  no  Reynolds  stresses  are  given. 

In  the  case  of  an  oscillating  airfoil  extensive  experiments  were  made  for  the 
determination  of  unsteady  airloads  by  Me  Croskey  and  his  collaborators.  A  summary 
and  discussion  of  these  experiments  are  given  in  C223.  Several  profile  shapes  were 
used  for  these  investigations  Unsteady  pressures  and  airloads  are  searched  for 
various  combinations  of  Mach  number  (up  to  .3)*  frequency  (k*=.  2)*  oscillation 
amplitude*  mean  incidence*  Reynolds  number  (up  to  4*0 00*000)  and  airfoil  shape 
(NACA  0012*  Vertol  vr-7,  NLR  7301).  No  velocities  or  detailed  turbulences  were 
measured  Four  different  flow  regimes  were  considered  :  no  stall*  stall  onset* 
light  and  deep  dynamic  stall. 

A  second  type  of  experiments  on  oscillating  airfoils  is  performed  by  Me  Alister  and 
L.  W.  Carr  C18* 193.  The  experiments  consist  in  visualisations  at  low  Reynolds 
number  (21*000)  of  the  dynamic  stalled  flow  around  a  pitching  NACA  0012  airfoil. 
The  airfoil  oscillated  from  0  to  20  degrees  incidence  et  k*.  29. 


A  last  type  of  experiment  which  can  be  mentioned  is  performed  by  CM.  Ho  and 
S.  H.  Chen  C12#133.  The  purpose  of  this  experiment  is  to  analyse  the  flow  in  the 
vicinity  of  the  trailing  edge  of  a  plunging  airfoil  in  a  steady  free  stream. 
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Measurements  were  made  with  5  X-wires  in  the  wake  of  the  airfoil. 

Hence*  in  the  case  of  a  pitching  airfoil  little  or  no  detailed  information  about 
velocity  and  turbulence  is  available  at  the  present  time.  It  should  be  observed 
that  even  in  the  steady  state  case*  such  information  is  not  profuse.  Data  which 
are  closely  related  with  the  present  investigations  are  presented  by  J.  C.  Yu  £343 
In  this  reference*  velocity  and  three  Reynolds  stress  profiles  are  presented  at  16 
chordwise  positions  in  the  vicinity  of  the  trailing  edge  of  a  NACA  63-012  airfoil 
<  x/c  from  .9  to  1,1). 

Finally  it  should  be  mentioned  that  methods  for  the  theoretical  prediction  of  the 
af orementioned  flow  types  are  in  full  development.  Some  recent  works  can  be 
mentioned.  S. B.  R.  Kottapalli  and  Q.  A.  Pierce  £173  presented  an  analytical  study  of 
the  drag  on  an  oscillating  airfoil.  The  turbulence  input  was  based  on  the 
Cebeci-Smith  eddy  viscosity  model.  T.  Cebeci  presented  a  paper  £33  dealing  with 
the  calculation  of  boundary  layer  characteristic s  on  an  oscillating  airfoil* 
including  flow  reversal.  J.  Shamroth  £303  applied  a  time  dependent  compressible 
turbulent  Navier-Stokes  analysis  to  the  oscillating  flow  field  problem.  The  test 
Reynolds  number  was  2*600*000.  Turbulence  models  are  based  on  a  turbulence  energy 
equation.  Close  to  the  present  investigations  is  the  work  of  Baker  and  Yu  £13 
where  calculated  detailed  turbulence  and  velocity  profiles  are  presented  and 
compared  with  Yu's  data  £343  This  is  however  limited  to  the  steady  state  case. 

Although  probably  incomplete*  the  present  review  situates  well  the  present  work 
between  previous  and  actual  work  done  in  the  field  of  unsteady  aerodynami c s.  The 
present  detailed  data  are  still  rare  although  required  by  the  actual  development  of 
computat ional  techniques. 


3.  EXPERIMENTS  SEX  UE 


3. 1  Introduction 


The  complete  experimental  set  up  is  summarized  on  fig  1.  The  oscillating  airfoil 
with  the  driving  mechanism  are  represented  on  top  of  the  figure.  A  pitot  tube  is 
mounted  on  the  probe  support  close  to  the  hot  wire  for  calibration  purposes.  A 
step  motor  allows  for  the  pitchwise  translation  of  the  probe  which  rotates  around 
its  axis  through  a  second  step  motor.  A  single  49  degrees  slanted  hot  wire  is 
mounted  on  the  probe  support.  This  rotating  wire  allows  for  the  determination  of 
mean  flow  as  well  as  all  non-zero  Reynolds  stresses. 


Four  signals  are  controlled  from  outside  the  windtunnel*  a 
the  driving  mechanism*  the  hot  wire  signal  and  the  probe 
controls.  The  rotational  position  is  kept  constant  during 
therefore  not  recorded  continuously  but  written  manually 
of  each  probe  translation.  Other  signals  are  recorded. 


triggering  signal  from 
translation  and  rotation 
probe  translation  and  is 
on  the  tape  at  beginning 


Sampled  and  digitized  data  are  simultanuously  written  on  tape*  while  one  of  the 
recorded  signals  can  be  observed  on  a  scope  after  DAC  conversion.  The  obtained 
records  are  next  forwarded  to  the  computer  centre  for  processing.  Measurement  of 
new  and  processing  of  previous  experiments  occur  simul tenuously*  which  represents 
an  important  time  gain. 


3.  2  Oscillating  airfoil 

A  60  cm  chord  NACA  0012  airfoil  is  used  for  the  experimental  set  up.  This  airfoil 
is  mounted  in  a  2  by  1  meter  test  section  in  one  of  the  VUB  blowdown  windtunnels* 
where  speeds  of  17  m/s  are  reached  at  a  free  stream  turbulence  level  less  than  .  3X 

Kraemer's  rule  was  used  to  design  a  tripping  wire  which  is  fixed  at  10%  chord 
distance  from  the  leading  edge.  The  diameter  of  the  wire  is  d  «  2  mm  where 
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The  airfoil  is  excited  by  a  300  Watt  motor  reductor  which  allows  for  an  oscillating 
frequency  up  to  8  Hz  at  5  degrees  oscillation  amplitude.  Sinusoidal  pulsation  is 
obtained  using  a  crank-connecting  rod  mechanism*  where  a  1  meter  long  connecting 
rod  allows  for  a  motion  which  is  sinusoidal  within  a  few  percent  distorsion.  The 
connecting  rod  is  close  to  the  lower  wind  tunnel  wall.  Speed  is  kept  constant 
whithin  one  oscillation  period  by  regulation  of  the  motor  shaft  velocity  (0-5000 
rpm).  The  crank  axis  delivers  a  one  bit  trigger  signal  which  is  recorded 
simul tenuously  on  tape. 


3. 3  Probe  support 

The  rotating  hot  wire  is  mounted  on  a  support  fixed  in  the  wind  tunnel.  A  step 
motor  allows  for  a  very  slow  and  constant  speed  translation  of  the  hot  wire  probe 
support  in  the  pitchwise  direction  across  the  wake  area.  The  resolution  of  the 
positioning  is  less  than  .002  7.  chord  (.01  mm)  in  the  pitchwise  direction  and  is 
about  .  1  7.  chord  in  the  chordwise  direction.  The  probe  itself  can  rotate  around 
its  axis  through  a  48  positions  step  motor  (only  12  are  used).  The  excentricity  of 
the  sensor  can  be  reduced  to  about  .01  X  chord  (.05  mm). 


3.  4  Hot  wire  apparatus 

The  used  hot  wire  is  a  45  degrees  slanted  DISA  gold  plated  hot  wire  type  55P02. 
The  sensor  diameter  is  .005  mm  and  the  length  is  1.25  mm.  The  wire  is  connected  to 
a  TSI  CTA  anemometer  bridge.  No  1  inear i sat i on  of  the  hot  wire  signal  is  used  in 
order  to  avoid  phase  distorsion. 

The  calibration  of  the  wire  occurs  in  situ  by  use  of  a  pitot  tube.  The  operating 
free  stream  velocity  when  calibrating  is  considered  to  have  a  constant  value  of 
11.  5  m/s*  although  small  variations  from  this  mean  value  are  observed  over  the  long 
measurement  periods.  This  approximation  is  done  in  order  to  have  a  coherent 
reference  value  over  all  the  experiments 


3  5  Data  acquisition  f ac i 1 i tu 

The  data  acquisition  facility  of  the  VUB  Fluid  Mechanics  Department  is  an  essential 
and  unique  facility  for  the  present  investigation. 

The  aim  of  this  facility  is  to  sample*  digitize  and  write  on  tape  up  to  8  analog  or 
digital  input  signals.  Each  analog  signal  (eg.  hot  wire  anemometer  voltage)  is 
filtered  and  DC  amplified  by  the  system  before  being  sampled*  digitized  through  a 
10-bit  ADC  and  recorded  on  tape.  All  acquisition  parameters  such  as  sample 
frequency*  cut-off  filters*  run  time*  formatting  of  digitized  data  are  simply  read 
in  from  a  tty*  which  makes  the  system  flexible  and  widely  applicable.  The  maximum 
data  transfer  rate  is  about  150  kHz  Comments  can  be  added  on  the  tape  from  the 
tty. 

A  constant  DC  offset  is  subtracted  from  the  hot  wire  anemometer  signal  before  being 
DC  amplified*  in  order  to  have  an  acceptable  resolution  for  the  fluctuating 
amp  1 i tudes 

A  1  bit  wing  position  reference  signal  is  combined  with  a  9  bit  probe  position 
coordinate  and  recorded  as  one  10  bit  signal  on  the  tape 

A  detailed  description  of  this  datalogger  can  be  found  in  C113 


4.  ROTATING  SI-ANTED  Hfll  «IR£  MEASUREMENT  TECHNIQUE 


4.  i  Intr9jte<;ix9.ii 

Two  methods  for  the  determinat ion  of  Reynolds  stresses  in  a  flow  field  are  mainly 
discerned.  In  the  first  an  X-wire  is  set  up  in  such  a  way  that  the  wires  are 
sensitive  to  the  selected  turbulent  quantities.  In  the  second  method  a  single 
normal  or  slanted  hot  wire  is  used  and  in  order  to  obtain  sensitivity  on  the 


6 


different  turbulence  velocities*  it  is  rotated  around  its  probe  support. 

The  latter  method  has  some  specific  advantages  over  the  X-wire  technique.  It 
allows  for  the  determination  of  froth  mean  flow  and  all  turbulent  quantities*  which 
is  not  the  case  with  an  X-wire.  Moreover*  a  less  complicated  electronic  set  up  is 
needed  and  the  data  are  related  to  a  much  smaller  area. 

The  method  as  used  in  the  present  study  has  its  base  in  the  works  of  Kool  and  De 
Grande  CIO* 163.  A  special  attention  was  nowever  paid  to  the  choice  of  the  probe 
positioning  in  order  to  avoid  prong  effects  as  much  as  possible  and  to  simulta- 
nuously  increase  the  resolution  on  the  different  quantities  to  be  measured. 

The  geometry  of  the  wire  is  defined  on  fig  2.  ip  and  ^  are  the  spherical 
coordinates  which  define  the  flow  direction  in  the  absolute  frame.  ao  denotes  the 
slanting  angle  and  c.p  defines  the  probe  rotation.  The  yaw  angle  0  is  defined  as 
the  angle  between  the  velocity  vector  and  the  plane  normal  to  the  wire.  This  angle 
can  be  expressed  as  function  of  a0»  ap*  *  and  if  through  the  relationship 

sin$  =  s-ino^cosw  cosw 

+  cosoc  cosw  sinw  cosoc  12) 

o  P 

+  cosoc  sinw  since 

o  p 


Eq  <2)  can  be  obtained  through  the  scalar  product  of  unit  vectors  chosen  along  the 
hot  wire  and  along  the  velocity  vector. 


4.  2  llffn.  flow  Pjtr.ome.fr  qr$ 

As  the  wire  is  rotated  around  its  support  axis  for  different  «p  values*  variations 
of  the  yaw  angle  and  hence  of  anemometer  voltage  are  observed.  From  these 
variations  it  is  possible  to  determine  the  flow  direction  and  magnitude. 

The  anemometer  voltage  E  and  the  'effective'  cooling  velocity  Geff  (<G  when  4 
non-zero)  are  related  through  King's  law  by 

E2  =  fi  +  B  oJff  (3) 

A*  B  and  C  are  determined  from  a  in  situ  calibration  using  a  pitot  tube.  The 
cooling  velocity  Geff  depends  on  the  velocity  magnitude  Q  and  the  yaw  angle  6  A 
simple  cosine  law  can  be  used  as  a  first  ap pr o x ima t i on  when  relating  Geff  to  G. 

Q  .  =  Q  cos<fr  (4) 

eft 


This  law  is  valid  in  absence  of  prong  effects  ftr  high  wire  length  over  diameter 
ratio's  (400)  and  for  values  of  less  than  about  45  degrees  C22.  At  higher  values 
of  f  (or  lower  1/d  ratio)*  the  forced  convection  no  more  obeys  this  law  and  a 
simple  correction  can  be  introduced  through  Bruun's  law  C21 

Qeff  -  Q  cos1"*  (5) 

This  relation  can  be  used  for  values  of  <*  up  to  70  degr  C23 

As  the  yaw  angle  goes  to  90  degr*  the  velocity  becomes  tangential  to  the  wire  and 
eqs  (4)  and  (5)  would  predict  Geff  *  0.  This  is  an  unrealistic  behaviour  since  the 
wire  is  still  cooled  by  the  flow  in  this  situation.  Since  such  high  values  of  yaw 
angle  are  expected  in  the  present  study*  eqs  (4)  and  O)  are  corrected  to  the 
classical  law 

2  2  21/2 

Qeff*  0  tco*  ♦  *  (61 

Kool  CIO* 161  introduced  corrections  to  eq  (6)  in  order  to  take  account  of  prong 
effects  and  asymmetry  in  the  wire  responses  These  effects  will  however  be  avoided 
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through  an  adequate  probe  positioning  and  by  continuous  calibration  for  all  probe 
positions. 

It  is  to  be  observed  that  eqs  (5)  and  (6)  are  closely  related  since 


_2  „2,  2 .  .  m 

Q  (cos 


fJ2.  .  .  .  2  m 

=  Q  (3  -  sin 

-  aQ  o 

, _  .2  ,  eft  „  Z . 

~  lQef*’«.=0  +  '  — “2„Vo  *  *ln  ♦ 

*  asm  4» 


(2) 


2  1/2 

Q  tl  -  m  sin 


Hence*  from  eqs  (6)  and  (7) 


k2= 


1  -  m 


(8) 


Eq  (7)  will  further  be  used  instead  of  eq  <6>»  for  conveniance 

From  the  measured  anemometer  voltages  E.  or  the  corresponding  values  of  Geff  at 
different  ap  values*  one  has  to  determine  the  best  values  of  **  ,  m  and  G  This 
is  done  by  expressing  the  least  squares  condition  as  follows 


2  2  2 

Z  C  (G^^/G)  +  m  sin  <p  -  13  is  minimum 

a 

P 


(9  i 


From  eq  (9)  it  follows  that 


l/Q2  2  </ff+  a,  2  sin2*  O 2ff  =  2  Q2ff 


(10) 


2  "*  2  4  2 

l/Q*~  Z  Q^^sin^tt*  +  it*  Z  sin  <p  -  Z  sin  <\ 


an 


Z  C(Q<?ff/Q)2+  m  sin2«t>  -  13  sir»4>  cos<|>  0 


(  12) 


2  C(Q  ff/Q)2+  m 

2  <3$ 

sin  $  -  1 3  s>in$  co*$»  ^  -  0 

f  1^) 

Since  eqs 
to  (13)  as 

(10)  and 
foil ows 

(11)  are  1  inear 

in  l/Q  and  m.  it  is  possible  to  express  eqs 

(10) 

f  ( $> ,  q. )  =  0 

(14) 

g  ($*  q< )  *  0 

(It.) 

Since  in  the  present  study  the  considered  mean  flow  can  be  considered  as  two 
dimensional#  at  least  one  of  the  two  angles  v  or  w  can  be  considered  as  known.  The 
surfaces  of  two  dimensional  flow  will  be  defined  by  c  “  v*  ,  where  v-  is  to  be 
chosen  (see  4.4).  Hence  eq  (12)  can  be  dropped  and  the  least  squares  fit  can  be 

reduced  to  the  solution  of  one  implicit  equation  f<v>*0.  This  is  a  periodic 

function  in  ^  which  in  general  has  4  roots  (sometimes  8)*  as  illustrated  on  fig  3 
Two  of  the  solutions  correspond  with  maximum  values  of  the  sum  of  squares  and  can 
be  rejected*  the  two  other  roots  lie  at  180  degrees  interval  and  hence  give  the 
flow  direction*  but  not  the  flow  sense,  which  was  to  be  expected  since  the  wire 

cannot  know  from  which  side  it  is  cooled  This  point  will  further  be  discussed 

section  7 
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A  first  approximation  for  y  when  solving  f(^)=0  can  be  found  from 


Z  (Q-Q)sinoc 


tg  y*  =  sin  - - - —  (16) 

Z  (Q-Q)cosoc 

P 

On  fig  4,  an  example  is  given  of  the  experimental  and  fitted  directional  responses 
of  Qeff  as  function  of  ,  figs  3  and  4  nre  based  on  the  same  test  data 


ion  of 


S.kre$.sfi.s 


Flow  turbulence  can  be  observed  by  a  hot  wire  thanks  to  its  high  frequency  response 
(up  to  20  kHz).  The  intensity  of  the  anemometer  voltage  fluctuations  is  quantified 
bg 

E '  =  £  (E-E)231/2  (17) 

where  E  denotes  the  time  averaged  value  of  E(t)  and  E'  the  rms  of  the  voltage 
fluctuations.  From  this  rms  value,  one  can  estimate  Qeff',  the  rms  of  the 
fluctuations  of  Qeff  Di f f er ent ia t i on  of  eq  (3)  yields 

^  ,  E2-  H  ,1/C  -1  2E  ^  (18) 
dQeff  =  (  — §—  }  CB  °E 


(dQeff)2  = 


,  E  -  R  ,2/C  -2 
B  ’ 


2E  2  2 

§  )2  (dE)2 


Neglecting  second  order  terms,  eq  (19)  may  be  rewritten  as 


Q'  -  < 

ef  f 


E2-  R  , 1/C  -1 


When  obtaining  eq  (19),  it  was  assumed  that  the  fluctuations  Qeff'  and  E'  are  small 
enough  compared  with  the  averaged  values  5el?  and  E.  This  linearisation  for 
'small'  turbulence  (<107. )  will  be  applied  in  all  the  following 

The  velocity  fluctuations  from  the  mean  velocity  Q  are  denoted  bg  u,  v  and  w>  where 
u  is  directed  in  the  mean  flow  direction,  while  v  and  w  are  defined  as  respectively 
tangential  and  normal  to  the  two  dimensional  plane  v  *  y'  (see  fig  2)  Hence, 
these  components  can  be  written  as,  according  to  the  linearisation 


v  -  Q  dy» 


w  =  Q  c  osy,  dy» 


Differentiation  of  eq  (6)  yields 


dQ^  ^  f  —  R^dQ  +  Q  (  —  dy«  +  —  dip  ) 


8R  ,  8R 

^  u  1  u 

R  U  +  rr -  V  +  -  X -  W 

u  3q>  cosy,  dip 


2  1/2 

R  «  (1  -  m  sin  t) 
u 


(30) 
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aa 

3y 


u  sinA  3* 

5  -  m  -p—  cos*  ^ 

u 


and 


dfi  .  . 

u  si  n*  .  3* 

=  _  m  cos*  5J 

u 


.  3* 

cos*  =  -  sina^sin*+‘  cosy>  -  cosa^sxruy  sin^  cosa^ 

+  cosa  cosa*  sina 
o  p 


3* 

cos*  ^  =  -  sxna^cos'V  sin>p  +  cosaQCOsy>  cos^  cosa^ 


From  eqs  (24)  to  (28),  one  can  write 


dQ  =  flu+flv+flw 
eff  u  v  w 


wh  ere 


0  2^,1/; 
fl  =  <  1  —  m  sin  *) 


(2t.) 


(2^) 


(22) 


(20) 


(29) 


( 30 ) 


fl 


v 


m 


sin* 

IT" 


u 


(cosa  cosa;  since  -  si  no  siny>  cosw 
o  ^  p  o  ^ 

-  coso-  siri'i-  sinw  cosa  ) 
o  +  p 


(31  ) 


B 


w 


m 


sin* 


R 


u 


(cosa  cosif  cosa  - 
o  p 


si  na^s  i  n>p) 


r 


) 


Squaring  and  time  averaging  of  eq  (29)  yields 


(dQ  )2  =  (fl  u)2  +  (fl  v)2  +  (fl  w)2 
eff  u  v  w 

+  2fl  fl  uv  +  2fl  R  uw  +  2fl  fl  vw 
u  v  u  w  v  w 


(33) 


According  to  the  1  inear i sa t i on  assumption#  this  equation  can  be  rewritten  as 


2  2  2  2  2  2  2 

Q'*  =  R  u'  +  fl v'*  +  R 

eff  u  v  w 

+  2R  R  uv  +  2R  fl  uw  +  2R  fl  vw 

U  V  U  W  V  w 


(34) 


In  eq  (34)  all  six  components  of  the  Reynolds  stress  tensor  appear#  while  all  the 
coefficients  and  the  1  h.  s  are  known  values  for  each  rotational  probe  position. 
Hence#  at  least  six  probe  positions  allow  for  the  determination  of  all  six  Reynolds 
stresses*  as  far  as  these  positions  are  chosen  adequately  In  practice#  more  than 
six  positions  are  used  and  the  best  values  of  the  stresses  are  obtained  by  a  least 
squares  fit.  This  fit  is  performed  as  follows. 


mi  jfcjfaiif 


v^. 
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Eq  (34)  is  rewritten  as 

2  6 

Q'Z*  -  S  a  s  =0  C35) 

•*t  j=l  ■>  J 

where  Sj  denote  tha  Reynolds  stresses  and  a j  the  corresponding  coefficients  in  eq 
(34).  The  least  squares  condition  implies 

'2  6  2 

2  ( Q  ^  -  2  a  s  )  is  mini mum  ( 36 ) 

•ff  ,  j  j 

ap  j  =  l 

This  is  satisfied  when 

2  V2ajV  °effJ  =  ° 

°P  J 


2  V2*jV  Qeff>  =  0 

GC 

P 


These  6  equations  can  be  written  as 


2  b .  s.  as  c  i  *  1  /  6 

j  u  j  j 


cb:cs>  -  <c> 


b .  =  a  a 
l  j  i  j 

*  a  * 
P 


S  -  2  aiQef  f 

a 

P 


The  stresses  Sj  are  finally  found  from 

CSJ  =  CB3-1<C>  <42) 

In  the  present  study  the  flow  is  basically  two  dimensional  and  uu  and  vw  can 
therefore  be  considered  as  neglectable.  Hence/  the  rank  of  CB2*  <S>  end  <C>  is 
reduced  to  4. 

On  fig  4  an  experimental  and  fitted  directional  response  of  Gef f '  is  illustrated 

4.  4  Oitllaai  B.ra&t  mttLti.ojiina 

When  choosing  the  probe  positioning  with  respect  to  the  expected  flow  directions/ 
three  parameters  are  to  be  determined  :  the  inclination  a  main  value  for  the 
flow  angle  ^  and  the  range  of  «->  over  which  the  probe  will  be  rotated.  Two 
criteria  are  to  be  considered  when  choosing  these  parameters  :  prong  effects  are 
to  be  avoided  as  much  as  possible  and  the  wire  must  be  sensitive  to  all  the 
searched  parameters. 

Prong  effects  can  best  be  avoided  by  choosing  the  probe  axis  parallel  to  the 
expected  velocity  direction  while  they  induce  important  errors  as  *  or  \l/ 
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reach  90  degrees  or  more  E2,  10/  162  (probe  axis  normal  to  velocity)  The  position 
is  however  a  singular  situation  since  both  determinants  of  eqs  <10>-<11)  for 
the  main  flow  and  of  EB3  for  the  Reynolds  stresses  become  zero.  Hence/  an 
intermediate  solution  is  to  be  chosen. 

This  choice  is  based  on  an  analysis  of  CB21.  Since  the  components  of  CB3  are  sums 
over  trigonometric  functions/  the  order  of  magnitude  of  its  components  is  1  Th^ 
ranges  of  $ ,  and  ap  will  be  chosen  in  such  a  way  that  al }  components  of  LB3~ 

have  this  same  order  of  magnitude  (10  different  components/  CB3  being  symmetric). 

An  analytical  expression  for  these  components  has  not  been  found  and  the  analysis 
is  therefore  based  on  plots  of  these  10  components  as  functions  of  f  and  O'.  On 
figs  5  to  9.  logaritmic  plots  of  the  moduli  of  the  components  are  made.  On  figs  5 

and  6/  one  can  observe  the  behaviour  of  CB2"1  as  v?  resp.  O'  are  increased  from  O  to 

45  degrees.  When  comparing  these  figures  it  appears  that  the  best  behaviour  is 
found  when  choosing  =*  O.  At  a  value  0  =  0  and  $  *  15  degr*  the  largest  component 
has  the  desired  order  of  magnitude  (fig  5)/  while  at  *  0  (fig  6)/  this  magnitude 
is  still  10  at  a  flow  angle  of  30  degrees.  At  higher  values  of  $  and  0#  it  is  no 

more  possible  to  rotate  the  wire  over  360  degrees  without  masking  the  wire  by  one 

of  the  prongs.  Probe  rotation  is  then  limited  to  180  degrees  in  order  to  avoid  too 
important  prong  effects  On  figs  7  and  8  are  plotted  the  CB3~^  components  at  180 
degr  rotation  for  a  probe  axis  perpendicular  to  the  flow  or  ^  =  90  degr).  From 
these  figures  it  appears  that  besides  the  prong  effcts.  a  bad  behaviour  is  found 
for  CB3~!  at  a  perpendicular  probe  set  up>  since  some  components  reach  values  of 
100  (fig  7)  and  1000  (fig  8) 

According  to  these  considerations,  it  was  decided  to  choose  f  as  close  to  zero  as 
is  allowed  by  the  magnitude  of  the  CBI-1  components.  This  value  is  chosen  as  t  = 
15  degrees.  On  fig  9/  it  is  shown  that  the  magnitude  of  CB3~!  remains  advantageous 
at  any  flow  angle/  except  around  &  =  42  degrees  (or  -42  degrees)*  where  some 
components  reach  values  of  10  and  100.  Hence  less  resolution  and  more  scatter  on 
some  results  can  be  expected  at  these  flow  angles 

The  largest  components  of  CB3'1  moreover  determine  which  Reynolds  stresses  will  be 
most  scattered  According  to  the  present  analysis*  it  was  found  that  most  scatter 
will  be  found  on  v'  and  next  on  w'. 


5  DATA  REDUCTION 


5.  1  Introduc  t  i  on 

A  main  problem  of  the  present  study  is  the  acquisition  and  processing  of  an 
enormous  amount  of  data.  About  109  bytes  have  been  recorded  (40  1600-bpi  2400  ft 
magnetic  tapes)  and  have  been  processed  on  a  170/750  CDC  computer. 

The  amount  of  data  to  be  recorded  was  chosen  as  a  compromise  between  minimum 
requirements  for  obtaining  an  acceptable  accuracy  on  the  results  at  one  hand/  and 
acceptable  measuring  times  and  computer  requirements  on  the  other  hand/  as  will  be 
outlined  in  the  following. 


5  2  Ensemb 1 e  average 

The  theory  described  in  section  4  is  valid  for  a  steady  flow  while  in  the  present 
i nvest igat i on  a  large  scale  fluctuation  on  E  is  impressed  by  the  periodic  motion  of 
the  blade.  Hence  averages  are  to  be  taken  at  times  when  main  values  of  £  and  E' 
are  identical,  which  can  be  performed  by  a  periodic  sampling  (ensemble  average) 
In  this  way  values  with  overbars  such  as  E  are  to  be  interpreted  as  periodic 
averaged  values  which  vary  periodically  with  time 

The  amount  of  periodic  samples  to  be  taken  is  determining  for  the  required 
measuring  times  and  is  therefore  reduced  to  an  absolute  minimum  of  30.  At  each  of 
these  30  osc i 1 lat ions,  about  250  samples  of 


El  -  A  (E-EO) 
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are  taken  at  a  constant  sample  frequency*  EO  being  a  constant  offset  voltage  and  A 
an  amplification  factor  A  and  EO  are  chosen  in  such  a  way  that  the  best 
resolution  is  obtained  in  the  data  acquisition  system  <8-bit)  It  is  to  be 
observed  that  the  hot  wire  is  calibrated  through  El  values  and  not  E*  in  order  to 
increase  the  ralibration  resolution  and  to  eliminate  A  and  EO  fluctuations 


3  3  Instant* ntous  profiles 

In  order  to  easily  reconstruct  instantaneous  wake  profiles*  the  probe  is  translated 
slowly  across  the  swept  wake  area  when  measuring  (see  fig  10).  The  translation 
speed  is  chosen  in  such  a  way  that  the  distance  covered  by  the  probe  during  30 
oscillations  is  small  enough  compared  with  variations  of  the  instantaneous  wake 
profile  within  this  area.  This  distance  is  13  7,  chord  or  about  1.  5  X  of  the 
observed  wake  thicknesses.  Hence  results  will  be  available  at  100/1.5  or  about  65 
points  within  the  instantaneous  wake  profiles,  which  allows  for  the  drawing  of 
quite  continuous  experimental  curves. 

Figure  10  shows  the  time  signals  as  they  are  recorded  on  digital  tape.  Ensemble 
average  over  30  oscillations  delivers  profiles  as  shown  on  fig  11  (with  much 
scatter  however).  A  cut  through  these  profiles  with  voltage/space  planes  delivers 
instantaneous  space  profiles  of  E(t)  and  E'(t)  Up  to  250  cuts  can  be  made*  but  in 
order  to  limit  the  computer  requirements*  only  46  are  taken  (each  7.  5  degr  phase 
shift). 

The  triggering  in  time  occurs  with  a  1  bit  reference  signal  delivered  by  the 
driving  mechanism  Triggering  occurs  with  an  accuracy  of  1/250  ot  .4  7,  of  the 
oscillation  period*  since  250  samples  of  this  signal  are  taken.  At  maximum  wing 
speed  this  corresponds  to  .An  7,  or  1.25  %  of  the  distance  swept  by  the  wing.  At 
maximum  oscillation  amplitude  <5  degr)*  this  corresponds  to  about  2  %  of  the  wake 
thickness*  which  is  close  to  the  1.5  7.  distance  over  which  the  data  are  taken. 
Hence*  instantaneous  profiles  with  about  65  significant  experimental  points  will  be 
available  at  48  times  during  one  oscillation  period. 


5  4  Conditional  averaging 

Averaging  over  30  events  is  insufficient  for  obtaining  results  without  important 
scatter.  More  accuracy  would  be  obtained  with  at  least  200  or  300  events*  but  this 
would  imply  about  10  times  more  measuring  time  (one  week  for  one  low  frequency 
experiment)  In  reality*  more  than  30  values  are  taken  into  account  when  searching 
for  the  turbulence  levels.  Since  much  data  are  available  in  space  as  well  as  in 
time*  it  is  possible  to  take  conditional  averages  over  small  space/tame  areas.  In 
reality  about  400  data  values  are  taken  into  account  in  this  way  in  order  to  obtain 
E ' .  Care  must  be  taken  however  in  order  to  not  detroy  significant  detailed 
information  and  the  used  averaging  procedures  which  are  summarized  in  section  55 
were  merely  obtained  by  trial  and  error 

The  most  helpful  step  when  averaging  is  illustrated  on  fig  12  and  is  performed  as 
follows.  The  48  instantaneous  profiles  are  translated  along  the  space  coordinate 
in  such  a  way  that  the  minimum  velocity  coordinates  coincide  for  all  the  profiles 
(see  fig  12).  When  a  time  cut  is  made  through  these  profiles*  one  obtains  a 
periodic  function  where  the  periodic  variations  are  rather  slow  compared  with  the 
scatter  which  can  therefore  easily  be  reduced  by  a  conditional  averaging*  without 
affecting  the  details  along  the  space  axis 


3.  3  SMfflffltrjI  (M  comoutT  raauiramant* 

The  different  steps  of  the  data  reduction  are  summarised  in  the  following. 

i)  Ensemble  average  from  magnetic  tape  data*  for  one  space  coordinate 

ii)  Weak  smoothing  of  E(t)  and  E'(t)  profiles 

iii)  48  phase  angles  are  selected  from  the  data  and  saved  on  mass  storage 

iv)  Repeated  for  each  space  coordin^t^ 


- - - -  —  — 
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v)  Repeated  for  each  probe  position  02) 

These  steps  require  about  1000  System  Seconds  on  the  170/750  CDC  computer  for  each 
experiment.  Data  are  reduced  from  50  Mbyte  to  about  200  10-bit  Kbyte. 


vi)  Instantaneous  E  and  E'  profiles  are  reconstructed  from  mass  storage  for  all 
probe  positions!  for  one  phase  angle. 

vii)  Weak  smoothing  of  these  profiles. 

viii)  Least  squares  fit  (see  4.3  and  4.4)  in  order  to  obtain  velocity#  flow  angle 
and  four  non-zero  Reynolds  stress  profiles.  These  profiles  are  defined  in  the 
absolute  frame. 

ix)  Results  are  saved  on  mass  storage 

x)  Repeated  for  each  of  the  48  phase  angles 

These  steps  require  about  110# 000  CM  words  and  200  System  Seconds  for  each 
exper iment. 


x i )  Profiles  are  read  from  mass  storage#  aligned  by  minimum  velocity  and 
information  which  lies  within  the  wake  is  windowed  (see  fig  12) 

xii)  The  obtained  profiles  are  smoothed  in  time  (fig  12)  eliminating  most  of  the 
scatter. 

xiii)  All  desired  output  quantities  are  derived  (see  section  7.1).  plotted  and/or 
saved  on  mass  storage  for  final  availability  Final  results  are  saved  on  magnetic 
tape. 


6.  SELECTION  OF  TEST  CASES 


The  selection  of  test  cases  was  mainly  based  on  previous  work  done  on  oscillating 
airfoils  C4> 12# 13.221  The  parameters  to  be  chosen  are  the  oscillation  amplitude# 
the  mean  airfoil  incidence  and  the  oscillation  frequency.  Chord  Reynolds  number  is 
kept  constant  at  about  300.000  All  test  are  done  at  .7#  2#  5#  10  and  20  X  chord 
downstream  of  the  trailing  edge.  All  done  experiments  are  listed  in  table  2 


6  l  QALLLLallaa  «mnlitud» 

Two  values  of  oscillation  amplitude  are  considered  The  maximum  value  used  is  5 
degrees#  this  value  has  been  used  in  the  past  in  oscillating  airfoil  experiments 
C222  At  this  amplitude  the  area  swept  by  the  wing  is  large  compared  with  the 
unstalled  wake  thickness  (13  X  chord  against  8  X)  and  the  incidence  varies  in  an 
important  way  (10  degrees).  Tests  were  also  made  at  a  low  amplitude  of  17  degrees 
(swept  distance  is  4.4  X  chord)  in  order  to  compare  results  at  similar  maximum  wing 
speeds 


6  2  Mean  incidence 

Considered  mean  incidences  are  0#  5#  9  and  10  degrees  Maximum  incidences  at 
maximum  oscillation  amplitude  are  resp  5#  10#  14  and  15  degrees  At  O  and  5  degr 
mean  incidence  no  stall  or  flow  reversal  is  expected.  At  9  degr  mean  incidence  (14 
degr  max  inc)  stall  onset  is  expected  according  to  the  measurements  of  Me  Crosfcey 
et  al  C223#  while  at  10  degr  mean  incidence  (15  degr  max  inc)  light  dynamic  stall 
is  present  Deeper  stall  is  not  considered  at  the  present  time 
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6.  3  Quillatjgp  ftiJutncu 

Three  situations  in  considered  steady  state#  'low'  and  'high'  frequency  'Low' 
frequency  means  that  the  passage  time  from  leading  to  trailing  edge  of  the  boundary 
layer  fluid  is  small  compared  with  the  oscillation  period  Hence#  this  situation 
can  be  considered  as  close  to  consecutive  steady  state  situations#  although 
significant  disturbances  due  to  the  wing  motion  are  found  At  'high'  frequency, 
the  passage  time  from  leading  to  trailing  edge  has  the  same  magnitude  as  the 
oscillation  period  and  large  phase  shifts  are  found  between  wing  oscillation  and 
turbulence  response  in  the  wake  and  in  the  rearmost  part  of  the  wing  boundary 
layers 

The  corresponding  frequency  coeffi^  jnts  are  36  and  1  to  1  23  Corresponding 
frequencies  2  16  and  6  1  to  7,4  Hz 

The  considered  frequencies  are  small  when  compared  with  the  'bursting'  frequency 
The  unsteady  effects  in  the  present  experiments  are  however  mainly  due  to  the 
strong  variation  in  incidence  and  the  propagation  of  the  disturbances  along  the 
airfoil#  whereas  in  experiments  at  fixed  incidence  no  similar  behaviour  is  present 
and  effects  can  best  related  to  something  as  the  'bursting  frequency'. 


7  RESULTS 


All  obtained  results  are  reported  on  figs  16  to  69  Continuous  curves  are  drawn 
since  about  100  experimental  points  are  available  for  one  curve,  allowing  for  the 
observation  of  small  details  The  results  are  drawn  each  45  degrees  phase  angle 
The  motion  of  the  blade  is  defined  through  the  instantaneous  incidence  i  as  follows 

cosw  f44) 

o  a 

where  i  *  instantaneous  incidence 

i  m  mean  incidence 
*  amp  1 1 1 ud e 
*  2*  ft  phase  angle 


*  *  0  always  corresponds  to  the  minimum  incidence  Data  are  shown  at  7  (  ,),  2 

(t').  5  <A>#  10  < f  1 )  and  20  <  )  V.  chord  downstream  of  the  trailing  edge 

The  data  as  shown  in  the  present  report  are  saved  on  magnetic  tape  and  are 
available  each  7  5  degrees  phase  angle 

Before  going  into  the  discussion  of  the  results#  details  will  be  given  of  the  exact 
definition  of  the  output  variables  and  the  way  used  to  calculate  them 


7  1  QkilUi  var labltt 

The  basic  data  as  detailed  in  section  4  are  defined  in  the  absolute  wind  tunnel 
reference  system  For  a  coherent  representation  of  the  results#  they  will  be 
transformed  to  a  relative  system  fixed  on  the  oscillating  wing 

The  'chordwiee'  direction  in  this  system  is  defined  as  the  i ns t antaneous  centerline 
of  tha  MAC A  profile  The  'pltchwise'  coordinate  is  normal  to  it  and  has  its  zero 
reference  on  it  (see  fig  14)  Th»  absolute  velocity  O  is  projected  in  this 
coordinate  system  through 
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V  as  Q  cos(w~i) 
c 


(4b> 


V^=r  Q  sin(qi— i)  (46) 


where  the  instantaneous  incidence  i  is  found  from 

eq  (44) 

(a  correction  due 

to 

non 

infinite  crank /connect ing 

rod  ratio 

is  however 

taken 

into  account).  The 

absolute 

and  relative  chordwise  velocities  are 

identical 

V 

(r)  ss  V 
c  c 

(42) 

while  the  relative  pitchwise 

vel oc i ty 

is  given  by 

( see 

14) 

V  (r)  = 

P 

. .  di 

Vp  ~  r  dt 

(48) 

a 

V  +  r  i  2-rrf 
P  a 

si  nip 

The  relative  flow  angle  i//(r) 

is  next 

found  from 

V  (r) 

tg 

c 

(49) 

Drawn  output  variables  are  Vc.  Vp ( r ) 

and  \if  <r  ) 

The  turbulence  quantities  u' 

.  v  ' »  w  ' 

and  uv  are 

also 

to 

be  transformed 

to 

the 

relative  i gitem  u Mr)  is  defined  as  the  relative  streamwise  fluctuating  velocity 

and  v'<r>  the  fluctuation  normal  to  it  The  spanwis*  fluctuation  w'  is  unchanged 

w '  ( r )  m  w '  ( t>0 ) 

u'(r)  and  v'(r)  can  be  obtained  as  follows 

u(r)  *  u  cos<$  -  v  sin4  (t>l) 


v(r)  *  u  *m6  ♦  v  cos<4  (o2) 

6  -  <p  -  if(r)  (t>3) 

Squaring  and  averaging  of  eqs  (51)  and  (52)  yields#  according  to  the  assumption  of 
1  l  near  nation 

,,  ,2  .2  2,  ,2  2  ,  ---  .  .  (t>4) 

u  (r)  *  u  cos  6  ♦  v'  sin  d  -  2uv  sin<*  cosd 


( r ) 


2  2  , 

sin  4 


2 

cos  4 


♦  2*jv  sin6  cosrf 


(t>5>) 


—  2  2  _  2  2  f  ^ ) 

uv(r)  ■  (u'  -  v  )  cosd  sin<4  ♦  uvic os  6  -  sin  4) 

Drawn  output  variablas  are  u'(r),  v'(r).  w '  and  uv(r)  It  is  to  be  observed  that 
far  downstream  A  is  small  and  no  differance  is  observed  between  relative  and 
absolute  values  of  the  Reynolds  stresses,  within  the  considered  oscillation  speeds 

It  is  possible  to  calculate  the  values  of  wake  displacement  and  momentum 
thichnesses  as  periodic  functions  from  the  detailed  data  These  values  are  defined 
as 

a’  •  ;**iv  v  i <*j 

-m  C  o*  C 


(32) 


-  16  - 


B  =  r  (V  -  V  )V  dy 
'  -CO  Coo  c  c 


h  =  /a  ' 

where  g  denotes  the  pitchwise  coordinate  In  theory/  Vc^  can  be  defined  without 
ambiguity/  but  in  practice  it  is  not  possible  to  perform  the  numeric  integration  in 
pitchwise  direction  far  enough  to  obtain  Vc  -*  Vc*  /  since  large  velocity  fluctua¬ 
tions  are  found  at  the  edge  of  the  viscous  wake  (see  fig  15).  In  order  to  obtain 
consistant  results  which  are  almost  independent  of  the  chosen  window  (yl#y2)i 
following  definition  was  introduced  (fig  13) 

y. 

6  =  S  tv  V  ) dy  +  f  (V  0-  V  ) dg  (60) 

y.  cl  c  y  c2  c 

1  m 

y_  y7 

&  =  /  (V  -  V  )V  dy  +  J  (V  V  )V  dy  (61) 

y,  cl  c  c  y  c  2  c  c 

1  m 

* 

Drawn  output  variables  are  6  i  0  and  H#  by  use  of  eqs  (59)/  (60)  and  (61). 

7.  2  »t«t»  rt»ul tt 

Steady  state  results  are  available  at  0/  3>  10/  and  13  degrees  incidence.  Detailed 
results  are  reported  on  figs  16  to  18  while  wake  thickness  parameters  are  shown  on 
figure  19.  Results  of  the  measurements  made  at  13  degrees  incidence  are  not 
reported  :  these  results  are  not  reliable  due  to  the  strong  turbulence  intensities 
(stalled  flow).  13  and  13  degrees  are  incidence  values  just  below  and  just  above 
the  static  stall  limit. 

The  results  as  they  are  reported  on  figs  16  to  18  were  obtained  without  smoothing/ 
except  the  v'.w'  and  uv  profiles  near  stall  (13  degr). 
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The  symmetric  results  at  0  degrees  incidence  can  be  compared  with  calculated  data 
of  Baker  et  al  C13  (fig  21)  On  this  figure  similar  double  peaks  in  u'  and  single 
peaks  in  v'  are  observed.  The  outer  parts  of  the  turbulent  profiles  are 
practically  unchanged  as  the  flow  goes  from  .7  to  20  X  chord  distance/  while  in  the 
central  part  a  decrease  in  all  turbulence  profiles  is  observed  over  a  width  which 
gradually  Increases  downstream. 

The  values  of  wake  thickness  and  shape  factor  of  the  static  test  cases  are  shown  on 
fig  19  and  are  also  listed  in  table  3.  Results  at  O  and  3  degrees  incidence  are 
quite  similar/  while  a  rapid  Increase  in  thicknesses  and  shape  factor  is  observed 
at  10  and  13  degrees  incidence/  in  particular  close  to  the  trailing  edge.  Momentum 
thicknesses  are  practically  constant  over  the  3  downtream  distances. 


7  3  Zero  flifjn  inc  idence  cases 

These  results  are  reported  on  figures  22  to  31, 

On  figure  22,  the  wake  thicknesses  are  first  shown  as  functions  of  time#  at  the  5 
considered  downtream  positions  At  low  frequency  <WL05)»  small  disturbances  are 
found  on  the  displacement  thickness.  Momentum  thikness  and  shape  factor  show 
larger  oscillations  with  a  maximum  shape  at  about  49  degr  phase  lead  to  maximum 
incidence  times  (0  and  180  degr  phase).  The  wake  thicknesses  are  very  close  to  the 
steady  state  values.  At  high  frequency  but  low  amplitude  (WH01>#  oscillations  are 
lightly  more  pronounced.  The  phase  lead  of  maximum  shape  is  increased  to  about  70 
degrees.  At  high  frequency  and  high  amplitude  (WH05)#  important  oscillations  and 
more  distorsion  are  found  in  the  thicknesses.  The  maxima  of  shape  factor  are 
increased  to  2.  close  to  the  trailing  edge. The  same  phase  lead  of  70  degrees  with 
respect  to  maximum  incidence  time  is  observed  at  low  and  high  oscillation  amplitude 
(WH01  and  WH05). 

Detailed  results  of  the  zero  incidence  test  cases  are  next  shown  on  figs  23  to  31 
These  results  are  shown  at  0#  45#  90#  135  and  180  degrees  phase  angle#  which 

corresponds  to  instantaneous  incidences  i  of  resp.  -5#  -3.53#  0#  3.53  and  5 

degrees.  Results  are  shown  at  increasing  incidence  <di>0>#  while  symmetric  results 
are  found  at  decreasing  incidence  Profiles  e c  O  and  100  degrees  phase  should  be 
completely  symmetric#  but  small  differences  are  however  found.  Since  these 
asymmetries  are  similar  in  all  exper  •  ,-r  at  zero  mean  incidence,  they  are 
believed  to  be  due  to  windtunnel  effe  t*  anti:  clearance  problems  in  the  driving 
mechanism. 

At  low  frequency  (figs  23  to  25#  WL05>#  results  are  in  general  close  to  the  steady 
state  results#  except  at  .  7  7,  downs^reatr  distance  where  the  wing  speed/flow  speed 
ratio  reaches  about  .25.  At  this  distance  details  in  u'  profiles  similar  to  those 
at  5  degrees  steady  state  incidence  are  observed  at  0  and  180  degrees  phase  angle# 
where  the  wing  speed  is  zero  (di=0>.  At  maximum  wing  speed  <90  degr  phase) 
detailed  differences  are  observed  Less  difference  is  observed  on  the  other 
profiles.  Central  peaks  are  found  again  in  the  v'  and  w'  profiles#  but  not  in  a 
reliable  way  since  they  are  not  present  at  0  degr  phase  angle. 

At  high  frequency  and  same  amplitude  (figs  26  to  28#  WH05)  differences  are 

increased  and  are  observed  up  to  10  7.  downstream  of  the  trailing  edge.  Velocity 
profiles  are  also  disturbed#  especially  at  pressure  side  close  to  the  trailing  edge 
(fig  26)  At  0  and  180  degr  phase  (di-0>#  detailed  turbulence  profiles  are  again 
rather  close  to  steady  state  as  well  as  low  frequency  profiles#  while  at  high 
speed#  important  differences  are  found  (compare  figs  24  and  27)  Central  peaks  in 
v#  and  w'  profiles  are  increased  at  maximum  wing  speed 

On  figs  29  to  31#  the  high  frequency  but  low  amplitude  test  case  is  shown.  Mo 
significant  detailed  differences  are  found#  although  the  maximum  wing  speed  at  90 
degr  phase  angle  is  the  same  as  in  the  low  frequency  and  high  amplitude  test  case 
(WL09)  Central  peaks  in  v'  and  w'  are  less  pronounced 


7  4  5  iPSS.POS  at 4H  amplitude  te.st  ca»es 
5  degrees  amplitude  cases  are  drawn  on  figs  32  to  50 

The  wake  thicknesses  are  first  reported  on  figure  32  The  maximum  shape  is  reached 
at  about  40#  70  and  90  degrees  phase  lead  to  maximum  incidence  time  at  frequency 
coefficients  of  resp  .36#  1.01  and  1.23  (resp  WL55#  WH55  and  WH51  >  The  maximum 
shape  is  almost  the  same  at  the  similar  maximum  wing  speeds  (WL55  and  WH51 )  while 
it  reaches  critical  values  at  high  wing  speed  (WH59#  H  *  2.6).  Hence  the  phase 
lead  depends  mainly  of  the  frequency  coefficient  while  the  maximum  shape  increases 
with  maximum  wing  speed 

The  most  important  oscillations  in  wake  thickness  are  found  at  low  frequency  (WL59) 
where  maximum  thickness  is  reached  at  about  60  degrees  phase  lag  after  maximum 
incidence  time  This  behaviour  was  to  be  expected#  according  to  the  history  of  the 
fluid  when  passing  from  leading  edge  to  trailing  edge#  since  this  fluid  has  seen 
the  maximum  incidence  some  time  before  reaching  the  region  just  behind  the  trailing 
edge  At  small  oscillation  amplitude  (WH5I)#  less  oscillations  are  observed#  which 


was  again  to  be  expected  since  the  absolute  level  of  incidence  does  not  vary  tou 
much  (3.3  to  6.7  against  0  to  10  degr). 

At  high  frequency  and  high  amplitude  (WH55)  a  quite  different  behaviour  of  the 
thicknesses  is  observed.  Oscillations  are  less  important  in  amplitude  (observe 
difference  in  scales)  and  maxima  of  h  are  shifted  to  the  increasing  incidence 
times  instead  of  the  decreasing  ones.  Maxima  of  shape  and  of  o*  are  almost  in 
phase#  while  at  low  frequency  they  were  shifted  over  180  degrees.  Maximum 
thicknesses  are  in  general  smaller  than  the  10  degrees  steady  results 

Detailed  profiles  are  next  shown  on  figs  33  to  50.  Detailed  results  are  drawn  at 
0#  45#  90  and  135  degrees  phase  angle  at  increasing  incidence  (di>0,  upper  figures) 
and  at  180,  225#  270  and  315  degr  phase  angle  at  decreasing  incidence  <di<0»  lower 
figures).  Corresponding  incidences  are  0#  1.47,  5.  8.53  and  10  degrees.  At  low 
frequency  (figs  33  to  38#  WL55)  and  zero  phase  angle  (zero  incidence#  zero  wing 
speed)#  profiles  are  found  which  are  rather  close  to  the  zero  incidence  steady 
state  ones.  At  maximum  wing  speed  however  (90  and  270  degr  phase)  an  important 
hysteresis  is  found  on  all  profiles#  which  is  in  accordance  with  the  behaviour  of 
the  wake  thicknesses  (fig  32)  Corr esp ond ing  to  the  instant  of  maximum  shape 
factor#  a  maximum  peak  in  the  u'  profiles  is  observed  (fig  35,  u'  up  to  147.)  close 
to  the  trailing  edge,  which  rapidly  decreases  downstream.  Mo  significant  peaks  are 
observed  on  the  other  profiles  Increased  values  of  uv  are  found  when  the  shape 
factor  is  maximum  at  2  and  5  X  downsteam  distance. 

At  high  frequency  and  high  amplitude  (figs  39  to  44#  WH55),  a  quite  different 
behaviour  if  found.  Profiles  at  zero  phase  angle  (zero  inc#  zero  speed)  are  rather 
close  to  the  5  degrees  steady  state  profiles#  which  can  be  explained  as  a  simple 
phase  shift  between  passage  along  the  airfoil  and  passage  at  trailing  edge  Flow 
reversal  is  approached  close  to  trailing  edge  when  shape  and  displacement  thickness 
are  maximum  (90  to  135  degr  phase).  High  corresponding  u'  and  w'  peaks  are 
observed  on  figs  41  and  44  These  peaks  close  to  trailing  edge  result  in  a  local 
smooth  increase  in  the  turbulence  profiles  at  20  7.  chord  distance  downstream  of 
trailing  edge  <u'  max  *  10  7.  at  20  7.  chord)  Mo  corresponding  peaks  are  found  in 
the  other  profiles  close  to  trailing  edge  but  from  5  to  20  X  downstream  distance# 
uv  is  increased  to  *.005  instead  of  -  002  to  -.003  elsewhere  (fig  42). 

The  highest  frequency  (WH51 >  but  low  amplitude  test  case  is  shown  on  figs  45  to  50. 
Some  detailed  differences  in  all  profiles  are  observed#  but  the  influence  of  the 
wing  oscillations  on  the  profile  details  is  in  general  small#  although  the 
oscillation  frequency  of  this  test  case  was  the  highest  one 


7  5  Sta 1 1  onset  test  cases 

Three  test  cases  were  considered  where  the  maximum  dynamic  incidence  is  near  the 
static  stall  limit  (14  degrees)  a  low  and  a  high  frequency  experiment  where  the 
maximum  incidence  is  15  degrees  (WI.A5#  WHA5) »  and  one  high  frequency  experiment  at 
14  degrees  maximum  incidence  (WH95> 

Wake  thicknesses  are  shown  on  fig  51  Although  the  stall  limit  is  lightly 

exceeded,  the  overall  qualitative  behaviour  is  similar  to  the  non-stalled  5  Degrees 
mean  incidence  cases.  Shape  factor  and  wake  thicknesses  are  however  increased  in 

an  important  way  The  maximum  shape  factor  in  the  high  frequency  case  is  again 

found  at  about  70  degrees  phase  lead  to  the  maximum  incidence  time#  while  at  low 

frequency#  these  maxima  are  spread  over  the  maximum  incidence  times,  or  when  the 

airfoil  is  above  the  static  stall  limit  Important  increases  in  shape  factor  are 

observed  at  high  frequency  up  to  3  6  for  WH95  and  4  1  for  WHA5 

On  figs  52  to  69  are  shown  the  detailed  results  of  the  test  cases  WLA5.  WH95  and 
WHA5.  Again#  a  quite  similar  qualitative  behaviour  as  in  the  5  degrees  mean 
incidence  test  cases  is  found  Peaks  and  turbulence  levels  are  however  increased 
At  low  frequency  (figs  52  to  57)»  flow  reversal  is  approached  around  30  degrees 
phase  lag  after  maximum  incidence  time  (fig  52  180  and  225  degrees  phase  angle), 

while  at  high  frequency  flow  reversal  is  touched  at  14  degrees  maximum  incidence 

(fig  58  :  WH95)  and  is  clearly  observed  at.  15  degrees  maximum  incidence  (fig  64 

WHA5).  This  flow  reversal  occurs  around  120  degrees  phase  angle#  at  increasing 

incidence.  The  flow  is  still  reverse  at  2X  downstream  distance#  while  at  5X  it 

has  disappeared  At  180  degrees  phase  angle#  when  the  Incidence  is  maximum  and 


above  the  static  stall  limit,  no  flow  reversal  is  observed 

High  turbulence  levels  are  observed  at  decreasing  incidence  in  the  low  frequency 
case  (figs  54  to  57)  and  these  profiles  can  be  considered  as  very  close  to  stall, 
according  to  fig  55  where  the  uv  profiles  are  shown  On  this  figure  a  behaviour  is 
found  which  is  believed  to  indicate  stall  onset  at  decreasing  incidence 

important  disturbances  are  observed  on  the  suction  side  of  the  uv  profiles  at  all 
downstream  distances.  uv  becomes  positive  at  pressure  side  of  the  wake  at  7  and  2 
V.  downstream  distance#  whereas  this  behaviour  disappears  further  downstream. 

At  high  frequency#  the  overall  turbulence  levels  are  more  uniform  over  the  whole 
oscillation  period.  Around  flow  reversal#  strong  up  and  downs  are  observed  in  the 
u'  profiles#  with  a  minimum  at  flow  reversal  (fig  66 )  The  corresponding  v'  and  w' 
profiles  show  a  different  behaviour  and  are  rather  increased  around  flow  reversal 
The  uv  profiles  (figs  61  and  67)  present  important  scatter  around  flow  reversal# 
probably  indicating  stall  onset. 

An  important  observation  here  is  that  even  at  low  frequency  or  when  flow  reversal 
is  present#  the  suction  side  of  the  wake  cannot  be  considered  as  completely 
stalled#  the  wake  thicknesses  and  turbulence  levels  being  small  when  compared  with 
steady  state  measurements  at  15  degrees  incidence.  Flow  reversal  is  detected  at 
high  wing  speeds  only 


8  CONCLUSIONS  AND  FUTURE  WORK 

In  the  present  report  the  technique  for  the  measurement  of  mean  flow  and  all 
non-zero  Reynolds  stresses  in  the  wake  of  an  oscillating  airfoil  is  detailed  All 
available  results  are  presented  and  discussed.  The  obtained  data  set  applies  on  a 
wide  range  of  oscillation  situations  at  several  positions  downstream  of  the 
trailing  edge  and  represents  an  extensive  contribution  for  direct  da ta-c ompar i s ons 

Following  qualitative  conclusions  from  the  above  discussion  can  moreover  be  made 

At  low  oscillation  amplitude  no  or  small  influence  is  observed  on  velocity  and 
turbulence  profiles#  within  the  considered  frequencies  and  in  the  non-stalled  case 

At  high  amplitude#  the  influence  of  the  oscillations  is  important.  A  phase  lead  in 
shape  factor  is  found  which  increases  with  the  oscillation  frequency  and  which  is 
the  same  for  all  mean  incidences,  except  at  stall  onset  The  maximum  thicknesses 
show  a  phase  lead  or  lag  to  the  maximum  incidence  time  depending  of  th"  mean 
incidence  and  frequency. 

Strong  disturbances  on  the  turbulence  are  present#  in  particular  at  high  wing 
speeds  close  to  the  trailing  edge  where  important  turbulence  peaks  are  observed 
which  result  in  an  overall  turbulence  increase  far  downstream. 

When  the  static  stall  limit  is  lightly  exceeded  (1  degree)  the  wake  cannot  be 
considered  as  stalled  at  frequencies  above  a  certain  frequency  which  is  less  than 
the  considered  low  frequency  <k-.3Q>.  Flow  reversal  is  detected  at  the  highest 
wing  speed  and  occurs  before  the  static  stall  limit  incidence  is  reached 

The  main  effects  of  the  oscillatory  motion  are  strong  phase  shifts  in  the 
turbulence  responses  and  detailed  disturbances-  in  particular  close  to  the  trailing 
edge 

In  the  near  future  the  present  investigations  will  be  extended  to  deeper  stall 
cases  and  the  measurements  will  be  extended  to  the  airfoil  boundary  layers  using  a 
sensor  mounted  on  the  airfoil  For  these  measurements  similar  test  conditions  will 
be  used  and  measurements  will  be  made  at  several  chordwise  positions  along  the 
oscillating  airfoil 
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